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I  .     INTRODUCTION 

A.  PURPOSE 

This  thesis  describes  the  modification,  completion  and 
initial  test  of  a  test  facility  designed  to  enable  the 
development  of  non- intrusive  measurement  techniques  within 
the  rotating  blade  row  of  a  transonic  turbine.   It  also 
details  the  computer  modeling  of  a  turbine  stator  flow  field 
for  validation  of  an  advanced  viscous  flow  code. 

B .  OVERVIEW 

The  performance  of  a  turbomachine  depends  heavily  on  the 
flow  fields  generated  within  the  machine  during  it's 
operation.   Measurement,  prediction  and  understanding  of 
conditions  in  the  flow  field  are  essential  to  optimized 
design  of  turbomachines .   More  efficiently  cooled  blading, 
enabled  by  advanced  computer  modeling  and  design,  can  result 
in  significant  improvements  in  service  life,  reliability  and 
lifecycle  cost. 

Installation  at  the  Naval  Postgraduate  School 
Turbopropulsion  Laboratory  (TPL)  of  the  cold- flow  test 
facility  repjrted  by  Studevan  [Ref .  1]  was  completed  as  part 
of  the  present  work.   The  facility  was  designed  to  enable 
tests  of  a  turbine  intended  to  power  the  High- Pressure  Fuel 
Turbopump  (HTJFTP)  (Fig.  1)  of  the  Space  Shuttle  Main  Engine 
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Figure  1.  Space  Shuttle  Main  Engine 


(SSME) .   The  turbine  was  designed  and  manufactured  by  Pratt 
&  Whitney  as  the  'Alternate  Turbopump  Development'  (ATD) 
model . 

The  purpose  of  the  new  Turbine  Test  Rig  (TTR) 
configuration  (Fig.  2)  at  TPL  is  to  facilitate  the 
advancement  of  non- intrusive  Laser -Doppler  Velocimetry  (LDV) 
techniques  required  to  validate  viscous  flow  simulation 
codes.   The  long  term  goal  of  this  research  is  to  use  LDV  to 
measure  the  flow  field  immediately  adjacent  to  the  endwall 
of  the  turbine  in  the  tip  clearance  region.   An  intermediate 
goal  is  to  determine  the  performance  of  the  single,  and 
ultimately  the  two -stage  turbine  operating  continuously  in 
the  TPL  facility.   Gaddis,  Hudson  and  Johnson  [Ref .  2] 
documented  the  exit  conditions  and  overall  performance  of 
the  ATD,  while  Hudson  and  others  [Ref.  3]  investigated  the 
original  turbine  model.   Both  of  these  tests  were  conducted 
in  a  short  deration  blow-down  facility  at  NASA  Marshall 
Space  Flight  Center . 

1.   Experiment 

Initial  measurements  of  the  ATD  will  determine  the 
power  output  of  the  single  stage  turbine.   The 
instrumentation  documented  herein  will  be  primarily  that 
concerned  with  the  safe  operation  and  control  of  the 
facility.   Provisions  have  been  made  for  initial  flow  field 
testing,  which  will  be  made  with  classical  intrusive 


Figure 


Turbine  Test  Rig 


techniques.   The  initial  tests  were  conducted  to  determine 
the  integrity  of  the  turbine  and  its  subsystems. 
2 .   Numerical  Simulation 

Numerical  simulations  of  the  flow  field  in  the  first 
stage  stator  of  the  SSME  HPFTP  turbine  were  accomplished 
using  two  codes  developed  by  Roderick  V.  Chima  called 
Turbomachinery  C-Grid  (TCGRID)  [Ref .  4],  and  Rotor  Viscous 
Code  3-D  (RVC3D)   [Ref.  5] .   TCGRID  was  used  for  the 
generation  of  the  3D  volume  grid  around  a  single  blade 
element  and  RVC3D  performed  the  numerical  solution  of  the 
f lowf ield . 


II.   FACILITIES 

A.   EXISTING  TEST  CONFIGURATION 

The  existing  configuration  of  the  TTR,  as  reported  by 
Studevan  [Ref .  1],  consisted  of: 

•  Inlet  Ducting 

•  ATD 

•  Shaft  and  Bearing  Housing 

•  Bearing  Housing  Support 

•  Dynamometer 

•  Bearing  Lubrication  System 

•  Controls  and  Instrumentation. 

All  of  the  above  listed  components  are  contained  within  the 
TTR  cell  with  the  exception  of  the  Controls  and 
Instrumentation  which  are  housed  in  a  control  room  adjacent 
to  the  cell.  The  interior  of  the  TTR  cell  may  be  seem  from 
the  control  room  through  an  explosionproof  window.  This 
allows  the  test  rig  to  be  operated  and  observed  from  a  safe 
environment.   The  test  arrangement  is  shown  in  Figure  2. 

During  operation,  high  pressure  air  flows  from  the  inlet 
ducting  through  the  inlet  strut  housing  and  into  the  ATD. 
In  the  ATD,  air  passes  through  the  annulus  of  the  first 
stage  stator  where  it  is  accelerated  and  turned  off  axis. 
The  air  then  impinges  on  the  first  stage  rotor  blades  which 


move  in  reaction  to  the  pressure,  turning  the  rotor  disk  and 
converting  the  energy  in  the  air  to  rotational  energy.   The 
air  then  exhausts  into  the  TTR  cell.   As  the  rotor  disk 
turns,  it  spins  the  main  shaft  in  the  bearing  housing.   The 
main  shaft,  via  the  quill  shaft,  is  coupled  to  the 
dynamometer  where  the  rotational  energy  is  absorbed. 

Shaft  speed  of  the  turbine  may  be  sensed  and  controlled 
by  the  dynamometer  system.   Configuration  of  the  speed 
control  system  is  unchanged  from  that  reported  by  Kane  [Ref . 
6] .   Speed  may  be  automatically  or  manually  controlled  by 
adjusting  the  water  flow  rate  through  the  dynamometer .   In 
the  automatic  mode,  a  magnetic  speed- sensor  signal  is 
compared  to  a  set  point,  and  a  controller  drives  the 
difference  to  zero.   The  dynamometer  is  instrumented  to 
provide  mechinical  power  output  via  torque  measurements. 
The  power  our.put  may  be  calculated  in  two  other  ways; 
namely,  usino  the  mass  flow  rate  and  temperature  rise  of  the 
water  through  the  dynamometer,  or  using  the  mass  flow  rate 
and  stagnation  temperature  drop  of  the  air  through  the 
turbine . 

B.   CONFIGURYTION  CHANGES 
1.   Bearings 

The  original  design  of  the  shaft  and  bearing  housing 
called  for  f^ur  high-precision  ball-bearing  units.   After 
coupling  the  shaft  and  bearing  housing  with  the  ATD,  the 


assembly  was  to  be  dynamically  balanced.   Balancing  proved 
to  be  impossible  and  the  assembly  was  taken  apart  for 
inspection.   Upon  disassembly,  one  of  the  bearing  units  was 
found  to  be  damaged.   In  order  to  simplify  the  system,  the 
outer  bearing  spacer  was  redesigned  so  that  only  two  ball 
bearing  units  were  necessary  (Fig.  3) .   Appendix  A  contains 
the  engineering  drawing. 

Duriig  reassembly  of  the  shaft  and  bearing  housing, 
the  shaft  was  cooled  in  a  refrigerator  to  62 *F  and  the 
bearings  were  slowly  heated  in  an  oil  bath  to  163 *F.   It  was 
noted  that  condensation  formed  on  the  shaft  when  it  was 
removed  from  the  refrigerator  and  it  was  wiped  down  prior  to 
inserting  the  bearing.   The  shaft  was  placed  on  end  on  a 
workbench.   When  the  bearing  reached  the  proper  temperature, 
it  was  removed  from  the  oil  bath  with  tongs  and  placed  on  a 
metal  plate.   From  the  plate,  it  was  immediately  picked  up 
in  gloved  hands  and  carefully  but  quickly  inserted  over  the 
shaft.   When  the  bearing  was  inserted  over  the  shaft  in  an 
attitude  level  with  the  workbench  it  did  not  bind  and  slid 
properly  into  position.   The  procedure  was  practiced  several 
times  with  the  damaged  bearing  before  an  attempt  was  made 
with  the  actual  hardware.   To  insert  the  shaft  and  bearing 
assembly  into  the  bearing  sleeve,  the  sleeve  was  preheated 
on  a  hot  plate  to  100 "F.   The  shaft  and  bearings  were  then 
pressed  into  the  bearing  sleeve  with  a  milling  machine  used 
as  a  press. 
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Figure  3.  Bearing  and  Shaft  Showing  Oil  Mist  Path 


2.  Oil-Mist  Cooling  System 

The  Portable  Turbine  Lubrication  Unit  described  by 
Studevan  [Ret.  1]  for  cooling  the  high-precision  ball 
bearings  was  replaced  by  an  oil -mist  lubrication  system. 
The  system  was  powered  by  shop  air  which  passed  through  a  5 
micron  filter  to  remove  particulate  matter  before  entering 
the  Norgren  Micro-Fog  Unit,  where  it  forced  a  small  stream 
of  atomized  oil  through  tubing  to  a  nozzle  (Fig.  4) .   The 
nozzle  was  inserted  into  a  hole  atop  the  downstream  end  of 
the  bearing  sleeve  (Fig.  3) .   The  oil/air  stream  flowed 
through  openings  in  the  outer  bearing  spacer  and  impinged  on 
the  ball  bearings.   It  then  passed  through  the  bearings  and 
out  two  ports  on  the  lower  upstream  and  downstream  ends  of 
the  bearing  sleeve.   This  was  an  open  system.   The  used  oil 
vented  to  waste  buckets  and  was  not  recirculated. 

3 .  Temperature  Monitoring  System 

Because  inadequate  cooling  could  lead  to  bearing 
degradation,  temperature- sensing  thermocouples  were  placed 
in  holes  through  the  bearing  housing  at  each  bearing  outer 
race.   These  were  wired  to  digital  output  meters  installed 
in  the  control  panel  in  the  TTR  control  room  and  were  used 
to  constantl/  monitor  the  bearing  temperature  and  insure 
that  the  cooling  system  was  operating  effectively.   The 
redesigned  bearing  housing  assembly  was  rebalanced 
successfully  at  approximately  2500  rpm  ,  with  the  oil -mist 
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Figure  4.  Oil  Mist  Cooling  System 
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lubrication  oystem  operating.   The  temperature-monitoring 
system  was  installed,  and  the  assembly  was  put  into  service. 
4 .   Vibration  Monitoring  System 

In  order  to  ensure  the  safety  of  operating  the  TTR, 
a  vibration- detection  system  was  installed  to  detect 
possible  out  -  of -balance  conditions  and  bearing  degradation. 
The  vibration- sensing  system  consisted  of  four 
accelerometers  and  a  signal  analyzer.   All  four 
accelerometers  were  attached  to  the  bearing-housing  support. 
One  was  upstream  at  the  top  of  the  housing  and  the  second 
was  upstream  but  displaced  90'  clockwise  looking  upstream. 
The  remaining  two  accelerometers  were  installed  on  the  two 
bear ing- housing  support  struts.   The  first  two  were  wired 
directly  into  a  Scientific -Atlanta  SD380  Signal  Analyzer  in 
the  control  coom.   The  two  accelerometers  installed  on  the 
support  struts  were  wired  into  the  control  panel,  and  to  an 
analog  Ballantine  Laboratories  RMS  indicator .   From  the 
analog  indicator,  they  could  be  connected  (if  needed)  to  the 
signal  analyzer  for  digital  processing  and  hard  copy  output. 

Prioi  to  powering  the  turbine  with  air  flow,  it  was 
analyzed  to  determine  natural  frequencies  so  that  resonance 
at  those  frequencies  could  be  avoided.   The  accelerometers 
were  installed  and  wired  to  the  signal  analyzer  before  any 
of  the  inlet  ducting  was  attached.   The  stationary  turbine 
and  bearing  housing  assembly  was  struck  with  a  modally  tuned 
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hammer  and  the  response  of  the  accelerometers  was  recorded 
by  the  signal  analyzer .   An  example  from  the  top 
accelerometex  may  be  seen  in  Figure  5 .   The  top  trace  shows 
the  real  time  response  and  the  bottom  shows  the  frequency 
spectrum  with  the  first  predominant  mode  at  67.5  Hz.   The 
assembly  was  then  turned  by  electric  motor  and  a  belt  drive 
at  varying  speeds  up  to  960  rpm  and  the  largest  amplitude 
responses  were  found  and  graphed,  as  shown  in  Figure  6 . 
Trace  A  is  from  the  top  accelerometer  and  shows  the  16  Hz 
speed  signal  along  with  its  higher  harmonics.   Trace  B  is 
from  the  accelerometer  at  90',  and  shows  similar 
characterist  Lcs . 

5.   Control  System 

The  control  system  in  the  TTR  control  room  included 
a  Fisher  -  Porter  electronic  set-point  controller  for  the 
outlet  valve  of  the  hydraulic  dynamometer .   This  was  used  as 
a  speed  regulator  with  manual  or  automatic  modes  and  was 
unchanged  from  the  configuration  described  by  Studevan  [Ref . 
1]  and  Kane  [Ref.  6] .   In  the  automatic  mode,  the  controller 
received  a  conditioned  feedback  signal  from  a  magnetic  speed 
sensor  in  ths  dynamometer .   The  controller  compared  the 
speed  signal  to  a  set-point  and  produced  a  corrective 
signal.   The  corrective  signal,  received  by  the  outlet 
valve,  controlled  the  water  level  in  the  dynamometer.   The 
speed  varied  in  response  to  the  resistance  produced  by  the 
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Figure  5.  Stationary  Shaft  Impact  Response 
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Figure  6.  Motor  Actuated  Rotational  Vibration  Response 
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water  and  the  energy  was  dissipated  into  the  water  as  heat. 
The  manual  inlet  control  valve  was  also  located  in  the  TTR 
control  room.   Manual  operation  of  the  speed  control  system 
was  used  during  all  the  tests  conducted  in  the  present 
study . 

6 .   Probe  Placement 

A  probe  access  hole  was  cut  through  the  outer  casing 
of  the  turbine  upstream  of  the  first -stage  stator .   The  hole 
location  corresponded  with  the  inlet  of  the  grid  designed  to 
be  used  in  the  computer  simulation  of  the  flow  field.   A 
probe  holder  was  designed  and  manufactured  to  hold  the  probe 
in  place  on  the  test  rig.   It  was  designed  to  be  used  in  any 
position  around  the  circumference  of  the  turbine.   However, 
a  new  access  hole  would  need  to  be  cut  in  order  to  use  it  in 
another  posicion.   An  engineering  drawing  of  the  probe 
holder  may  be  found  in  Appendix  A. 

Combination  probes  (numbers  2  and  3)  were  placed  in  the 
inlet  duct  through  holes  cut  in  the  piping.   They  were 
positioned  at  the  same  axial  station  to  measure  stagnation 
temperature  and  stagnation  pressure  into  the  turbine. 
Combination  probe  number  4  was  clamped  to  the  table 
immediately  downstream  of  the  turbine  to  measure  exit 
stagnation  temperature  and  pressure. 
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III.   EXPERIMENT 

A.   EXPERIMENTAL  PROCEDURE 

Before  any  measurements  were  attempted  in  the  TTR,  the 
shroud  inserts  between  the  blade  tips  and  the  turbine  outer 
casing  were  removed.   With  the  shroud  inserts  removed,  there 
was  no  danger  of  blade  tips  rubbing  due  to  misalignment. 
Initial  system  validation  runs  documented  herein  were 
conducted  with  the  shroud  inserts  removed.   The  TTR  was 
operated  twije  to  gain  experience,  test  safe  operating 
methods  and   shake-down'  the  facility.   During  the  first 
shake-down  run,  speed  control  and  vibration  response  were 
tested  under  various  air  and  water  flow  settings.   It  was 
noted  that  if  water  flow  was  sufficiently  reduced  by  closing 
the  outlet  valve  to  the  dynamometer,  the  speed  of  the 
turbine  climbed  rapidly.   On  inspection  of  the  TTR  after  the 
run,  the  dynamometer  was  found  to  be  very  hot.   During 
inspection  of  the  water -circulation  system,  water  was  found 
to  have  overflowed  onto  the  floor  of  the  pumping  station  due 
to  surging  developed  by  changing  demand  to  the  dynamometer . 
It  was  decided  that  the  water  flow  through  the  dynamometer 
should  be  set  at  the  start  of  the  run  to  keep  the 
dynamometer  cool  and  eliminate  surging  of  the  water 
circulation  system. 
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For  the  second  shake-down  run,  a  specific  procedure  was 
adopted  to  set  the  water  and  air  flow  at  the  start  of  the 
run.   In  order  to  keep  sufficient  and  consistent  water  flow 
through  the  dynamometer,  the  inlet  valve  to  the  dynamometer 
was  set  to  80%  open  and  the  outlet  valve  was  set  to  50% 
open.   For  mass  flow  to  the  TTR  to  be  properly  calculated, 
all  the  air  flowing  through  the  metering  orifice  had  to  go 
to  the  turbine  (Fig.  7) .   To  ensure  that  all  air  flow  went 
to  the  turbine,  the  number -two  dump  valve  had  to  be 
completely  closed.   After  the  TTR  inlet  valve  was  set  to  20% 
open,  the  number -two  dump  was  closed,  then  the  number -one 
dump  valve  was  set  to  40%  open.   Further  manipulation  of  the 
mass  flow  to  change  the  speed  of  the  turbine  involved 
opening  the  TTR  inlet  valve  to  increase  speed  or  opening  the 
number -one  dump  to  decrease  speed.   It  was  necessary  to  have 
some  air  venting  through  dump  valve  number  one  in  order  to 
keep  the  total  flow  rate  sufficiently  high  or  the  Allis- 
Chalmers  compressor  (Fig.  7)  could  surge  due  to  excessive 
back  pressure. 

During  the  experimental  run,  the  procedure  described 
above  for  setting  water  and  air  flow  was  followed.   Pauses 
were  made  during  the  acceleration  and  system  settings  were 
noted.   At  eich  pause,  the  speed  was  allowed  to  stabilize 
and  measurements  were  made  of  temperatures  and  pressures  for 
the  power  output  and  mass  flow  calculations.   Measurements 
relevant  to  power  calculations  are  shown  in  Table  I. 
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Figure  7.  Compressed  Air  System 
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TABLE  I .  ABBREVIATED  TTR  DATA 


Speed 
(RPM) 

(*F) 

(in.H20) 

Pi 

(in.  Hg) 

TT4 
CF) 

TT3 
CF) 

880 

81 

0 

37  .35 

77 

81 

1420 

88 

18  .438 

38  .57 

84 

88 

1560 

92 

19 .000 

38.98 

86 

92 

2100 

94 

16  .625 

41  .42 

87 

94 

3700 

97 

24  .500 

53  .23 

84 

97 

5950 

105 

33  .750 

49  .56 

76 

105 

7900 

111 

17  .  375 

55  .06 

71 

111 

B.   MASS  FLOW  AND  POWER  CALCULATIONS 

The  turbine  air  mass -flow  rate  was  metered  by  an  orifice 
plate  (Fig.  7)  and  was  calculated  using  the  method  described 
by  Vavra  [Ref .  16]  and  summarized  in  Appendix  B.   The  power 
output  of  the  turbine  was  calculated  using  the  stagnation 
temperature  drop  of  the  air  across  the  turbine,  and  the  mass 
flow  rate.   Table  II  shows  calculated  values  of  mass -flow 
rate,  specific  power  out  of  the  turbine,  and  efficiency. 
Details  of  these  calculations  may  be  found  in  Appendix  C  and 
a  complete  table  of  measurements  is  contained  in  Appendix  D. 
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TABLE  II.  TURBINE  PERFORMANCE  (SHROUDS  REMOVED) 


Speed 
(RPM) 

Mass  Flow 

Rate 

(lbm/s) 

Specific 

Power 
(Hp/lbm) 

Efficiency 
(%) 

880 

- 

- 

14  .96 

1420 

7  .51 

1.35 

7  .47 

1560 

7  .63 

2  .03 

9  .98 

2100 

7  .37 

2  .36 

6  .48 

3700 

10  .  08 

4  .40 

3  .86 

5950 

11  .27 

9  .  81 

3  .49 

7900 

8  .57 

13  .53 

3  .23 
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IV.   NUMERICAL  SIMULATION 

A.   GRID  GENERATION 

Numerical  simulation  of  a  flow  field  is  highly  dependent 
on  the  grid  used  to  define  the  geometry,  and  its  resolution. 
Many  routines  have  been  developed  to  generate  grids  within 
and  about  realistic  surfaces  that  bound  fluid  flows.   The 
computer  program  used  to  develop  the  grids  used  in  this 
investigation,  Turbomachinery  C-Grid  (TCGRID) ,  generates 
three-dimensional  C-  or  H-type  volume  grids  for 
turbomachinery.   Grid  files  are  compatible  with  RVC3D  and 
PLOT3D  [Ref .  7]  which  were  used  extensively  in  the  flow 
solution  and  visualization. 

TCGRID  generates  a  three-dimensional  computational 

volume  grid  in  the  following  manner: 

1.  A  course,  equally  -  spaced  meridional  grid  is  generated 
between  the  supplied  hub  and  tip. 

2  .  The  blade  geometry  is  interpolated  onto  the  meridional 
grid . 

3.  2-D  blade-to-blade  grids  are  generated  along  the 
meridional  grid  lines  in  (m, rbar* theta)  coordinates 
using  an  old  version  of  the  Sorenson  GRAPE  code 

[Ref .  3] . 

4.  The  (m, rbar * theta)  coordinates  are  transformed  to 

(z,r, theta)  coordinates. 

5.  The  2-D  grids  are  reclustered  spanwise  to  make  a  full 
3-D  grid. 
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6.  Finally,  the  (z,r,theta)  coordinates  are  transformed 
back  t;  (x,y,z)  and  written  in  PL0T3D  format. 


Hub,  tip  and  blade  geometries  must  be  in  MERIDL  [Ref .  9] 
format . 

The  hot -flow  blade  geometry  data  were  received  from  NASA 
Lewis  Research  Center  [Ref.  10]  in  (z,r,theta)  coordinates 
appropriate  for  use  in  TCGRID .   Since  the  Naval  Postgraduate 
School  facility  was  cold-flow,  a  99%  thermal  -  shrink  factor 
had  to  be  applied  to  all  linear  dimensions.   The  FORTRAN 
code  written  to  accomplish  this  may  be  seen  in  Appendix  D. 

Grid  generation  and  sizing  is  a  subjective  exercise  in 
which  there  is  no  single  'right'  answer.   Often  the  success 
of  a  grid  remains  unknown  until  it  is  used  with  a  flow 
solver.  In  this  investigation,  over  one  hundred  flow 
solutions  were  attempted  with  various  grid  sizes  before  a 
successful  combination  was  found.   It  was  understood  that 
the  investigation  required  a  fully  three-dimensional  viscous 
flow  solution.   Because  of  the  resolution  needed  in  the 
boundary  layer,  the  grid  had  to  be  very  fine.   However  ,\  the 
successful  grid  contained  more  than  twice  the  number  of  grid 
points  in  the  spanwise  direction  than  originally  thought. 
The  final  grid  size  used  for  the  flow  solutions  was  151  x  31 
x  45.   Chima  [Ref.  11]  recommended  15  grid  points  in  each 
endwall  boundary  layer  for  good  resolution.   An  example  of 
the  TCGRID  input  file  may  be  found  in  Appendix  D. 
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Figures  8  through  11  show  C- grids  developed  by  TCGRID  for 
the  hub,  midspan  and  tip. 

This  geometry  was  difficult  to  model  because  of  the 
endwall  contour  (Fig.  8) .   The  hub  and  tip  geometry  was 
defined  in  tne  input  file  with  (z,r)  coordinates.   Then  the 
stator  contour  was  defined  in  five  sets  of  (z,theta,r) 
coordinates.   Each  set  of  coordinates  was  at  a  constant 
radius.   Since  the  endwall  radius  reduced  at  the  stator  tip, 
some  of  the  coordinates  were  defined  outside  the  endwall. 
TCGRID  calculates  the  intersection  of  the  stator  coordinates 
with  the  endwall  and  generates  the  grid  defined  by  the 
intersection  contour .   Figures  9  through  11  show  grids 
developed  for  hub,  mean  and  tip  sections.   Figure  12  shows  a 
three-dimensional  view  of  the  final  grid.   The  curved 
trailing  edge  shown  in  this  figure  also  complicated  the  grid 
generation  procedure. 

B.   FLOW  SOLUTION 
1 .   RVC3D 

Once  a  proposed  grid  has  been  developed,  it  can  be 
used  by  a  flow  solver  to  simulate  a  flow.   RVC3D  (Rotor 
Viscous  Code  3-D)  is  a  computer  code  developed  at  NASA  Lewis 
Research  Center  for  the  analysis  of  fully  three-dimensional 
viscous  flows  in  turbomachinery .   The  code  solves  the  thin- 
layer  Navier -Stokes  equations  with  an  explicit  finite- 
difference  technique.   It  is  applicable  to  annular  blade 
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Figure  8.  Hub  and  Endwall  Contour 
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Figure  9.  C-Grid  at  Stator  Hub 
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Figure  10.  C-Grid  at  Stator  Midspan 
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Figure  11.  C-Grid  at  Stator  Tip 


28 


Figure  12.  Three-Dimensional  Stator  Grid 
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rows  or  linear  cascades  and  includes  two  turbulence  models 
along  with  a  simple  tip  clearance  model.   Mathematical 
formulation  of  RVC3D  is  described  in  references  12  and  13 . 

The  code  solves  the  Navier -Stokes  equations  in 
Cartesian  coordinates  with  rotation  about  the  x-axis.   The 
equations  are  mapped  to  a  general  body -fitted  coordinate 
system  (Fig.  13) .   Turbulence  is  modelled  using  a  three- 
dimensional  adaptation  of  the  Baldwin- Lomax  model  [Ref .  14] 
or  the  Cebeci -Smith  model  [Ref.  15] .   A  sample  RVC3D  input 
file  may  be  seen  in  Appendix  D. 
2 .   Solution  Method 

Simulations  were  run  on  two  different  computer 
systems;  a  Silicon  Graphics  Iris  workstation  and  a  Cray  Y-MP 
EL-98  supercomputer.   Both  systems  reside  at  the  Naval 
Postgraduate  School.   Solutions  were  obtained  for  two  inlet 
velocity  profiles  on  each  machine.   The  same  version  of 
RVC3D  was  loaded  onto  the  two  computers  and  identical  input 
files  were  used. 

First,  using  the  two  computers  for  time  efficiency, 
an  inlet  velocity  profile  to  simulate  a  fully  developed 
turbulent  flow  was  found  using  the  Iris  while,  concurrently, 
a  profile  simulating  a  slug  flow  inviscid  profile  was  found 
using  the  Cray.   The  profiles  were  developed  by  modifying 
the  RVC3D  input  files  and  running  100  iterations  at  a  time 
to  find  the  inlet  velocity  profile.   After  each  run,  the 
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Figure  13.  General  Body-Fitted  Coordinate  System 
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input  file  was  modified  and  run  again.   The  process  was 
repeated  until  a  satisfactory  profile  was  obtained.   It  was 
assumed  that  the  two  computers  would  yield  identical 
results . 

After  several  runs  on  each  computer,  it  was  noted 
that  when  similar  values  were  used  in  the  input  files  for 
the  two  machines,  the  velocity  profiles  that  were  calculated 
were  vastly  different.   The  turbulent  velocity  profile  was 
run  on  the  Iris  until  the  solution  converged  and  the  same 
was  done  on  the  Cray  with  the  slug  profile.   Then  the  input 
files  were  exchanged  and  the  simulations  were  run  again. 
All  four  solutions  were  plotted  for  comparison. 
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V.   RESULTS  AND  DISCUSSION 

A.   EXPERIMENT 

The  data  shown  in  Appendix  C  yield  curves  shown  in 
Figures  14  and  15.   The  curves  confirm  general  trends  that 
were  expected  from  the  experiment.   The  data  were  taken 
without  the  intention  of  producing  a  turbine  map  but  to 
validate  the  data  acquisition  methodology  and  to  check  the 
integrity  of  the  turbine  and  its  subsystems.   The  data  show 
that  the  syscem  will  work  but  that  it  is  inconsistent.   It 
was  noted  that  the  water  manometers  used  to  measure  pressure 
drop  across  the  metering  orifice  were  contaminated.   This 
explains  the  inconsistent  values  noted  in  Appendix  C  as  P31, 
P32,  P33  and  P34 . 

The  difference  in  the  stagnation  temperature  readings 
indicate  that  the  flow  may  not  have  been  circumf erentially 
uniform.   The  stagnation  pressure  values  PT2  are  obviously 
erroneous  and  the  data  obtained  with  that  combination  probe 
were  not  used. 

The  low  efficiencies  shown  in  Figure  15  were  expected. 
Since  the  shroud  inserts  were  removed  for  these  runs,  it  was 
expected  that  much  of  the  air  flow  would  pass  between  the 
blade  tips  and  the  outer  case.   The  air  passing  around  the 
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Figure  14.  Specific  Power  and  Mass  Flow  vs  Speed 
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Figure  15.  Efficiency  vs  Speed 
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rotor  blades  did  not  contribute  to  the  power  output  of  the 
turbine  and  significantly  reduced  the  efficiency. 

Figure  16  shows  the  speed  sensor  response  of  the  TTR  at 
800  rpm.   The  number  in  the  lower  left  hand  corner,  measured 
in  Hertz,  comes  from  the  rotational  speed.   But,  since  the 
magnetic  speed  sensor  on  the  dynamometer  detects  the  passage 
of  teeth  on  a  30  tooth  gear,  the  rpm  of  the  system  is  twice 
the  frequency  shown.   Figure  17  shows  the  speed  sensor 
response  at  7925  rpm. 

B.   COMPUTATIONAL  FLUID  DYNAMICS 

As  noted  previously,  running  RVC3D  with  the  same  input 
files  on  the  two  computers  resulted  in  solutions  with 
clearly  different  inlet  velocity  profiles.   Figure  18  shows 
the  two  turbulent  inlet  velocity  profiles  from  one  input 
file,  and  Figure  19  shows  the  slug  flow  profiles  from  the 
other.   The  differences  in  inlet  profiles  generated  in  the 
course  of  the  solution,  show  that  the  inlet  velocity  profile 
is  sensitive  to  either  the  architecture  or  the  precision  of 
the  particular  computer. 

Figures  20  and  21  show  the  differences  in  inlet  velocity 
profiles  using  the  two  input  files  on  the  same  computer . 
The  similarity  of  the  graphs  shows  that  the  output  inlet 
velocity  profile  is  not  very  sensitive  to  changes  in  input 
velocity  profile.   Comparison  of  Figures  18  and  19  with 
Figures  20  and  21  indicates  that  the  output  inlet  velocity 
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Figure  16.  Low  Speed  Vibration  Response 
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profile  is  more  sensitive  to  computer  than  to  input  velocity- 
profile  . 

Figures  22  and  23  show  the  midspan  chordwise  variation 
in  normalized  pressure  distribution  of  solutions  with  the 
same  input  velocity  profiles  on  different  computers. 
Figures  24  and  25  show  the  same  results  with  the  different 
inlet  profiles  on  the  same  machine.   Comparison  of  these 
four  graphs  shows  that  all  the  variations  gave  identical 
pressure  distributions.   Therefore,  the  pressure 
distribution  is  sensitive  to  neither  inlet  velocity  profile 
nor  computer . 

Figures  26  through  29  are  PL0T3D  plots  of  midspan  Mach 
number  contours.   The  first  two  show  the  Cray  solutions  of 
the  turbulent  and  slug  flow  inlet  velocity  profiles.   The 
latter  two  show  Iris  workstation  solutions  for  the  same 
input  files.   Again,  in  spite  of  the  differences  in  the 
input  files  and  the  velocity  profiles,  all  four  solutions 
are  identical .   The  plots  show  the  flow  accelerating  on  both 
the  suction  and  pressure  surfaces  of  the  stator  with  maximum 
Mach  number  in  the  area  of  the  throat  as  expected.   Boundary 
layer  growth  as  well  as  wake  shape  and  orientation  are 
evident  in  the  plots. 

It  is  important  to  note  that  the  wake,  in  yellow,  does 
not  follow  the  contour  of  the  grid  but  follows  the  exit 
angle  of  the  stator.   The  cut  that  joins  the  pressure  and 
suction  sides  of  the  grid  may  be  seen  in  Figure  10.   The 
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Figure  22.  Normalized  Pressure  -  Turbulent 
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Figure  23.  Normalized  Pressure  -  Slug  Flow 
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Figure  24.  Normalized  Pressure  -  Cray 
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Figure  25.  Normalized  Pressure  -  Iris 


4  7 


Figure  26.  Mach  Contours  -  Cray  -  Turbulent 
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Figure  27 .  Mach  Contours  -  Cray  -  Slug  Flow 
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Figure  28.  Mach  Contours  -  Iris  -  Turbulent 
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Figure  29.  Mach  Contours  -  Iris  -  Slug  Flow 
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fact  that  the  wake  does  not  follow  this  cut  indicates  that 
the  solution  is  somewhat  independent  of  grid  geometry  and 
that  continu'.ty  is  satisfied  at  each  point  in  the  flow. 

Figures  ."0  through  33  the  total  pressure  contours  across 
the  exit  pla  le  of  the  stator .   These  plots  are  also  nearly 
identical  but  small  variations  do  exist  in  the  upper  left 
corner  of  the  exit  plane.   The  turbulent  solutions,  Figures 
30  and  32 ,  are  different  in  this  region  from  the  slug  flow 
solutions,  Figures  31  and  33.   This  would  indicate  that  the 
solution  is  somewhat  dependent  on  inlet  velocity  profile. 

In  these  figures,  it  is  also  possible  to  see  the 
variation  in  wake  location  along  the  span.   The  yellow  area 
at  the  edges  of  the  plots  shows  that  the  wake  bows  toward 
the  pressure  side  at  midspan.   This  may  be  explained  by  the 
bowed  shape  of  the  trailing  edge  of  the  stator .   As  the  flow 
turns  around  the  stator,  the  extended  arclength  along  the 
surface  in  the  region  of  the  midspan  may  turn  the  flow 
farther  than  the  shorter  arclength  along  the  stator  at  the 
hub  and  tip. 

In  the  boundary  layer  at  the  hub  there  is  a  small  region 
of  blue  and  /ellow  in  the  area  of  the  grid  cut.   Although 
not  a  discon  linuity ,  this  locus  coincident  with  the  grid  cut 
may  indicate  some  solution  dependence  on  grid  geometry. 

An  example  of  the  residual  histories  of  two  solutions  is 
shown  in  Figure  34 .   The  plot  shows  that  the  residuals 
matched  exactly  for  the  first  3500  iterations.   The  Iris 
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Figure  30.  Total  Pressure  -  Cray  -  Turbulent 
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Figure  31.  Total  Pressure  -  Cray  -  Slug  Flow 
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Figure  32.  Total  Pressure  -  Iris  -  Turbulent 


Figure  33.  Total  Pressure  -  Iris  -  Slug  Flow 
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Figure  34.  Residual  History  Example 
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solution  was  continued  beyond  3500  iterations  in  order  to 
see  the  effect  of  varying  the  inlet  velocity  profile.   Each 
of  the  variations  of  the  input  file  is  marked  by  a 
significant  rise  in  residuals.   However,  after  approximately 
500  iterations  each  solution  converged  to  nearly  the  same 
residual  value.   Even  after  seven  changes  in  the  input  file 
to  vary  the  inlet  velocity  profile,  the  residuals  converged 
to  four  orders  of  magnitude  below  their  start.   This  shows 
that  the  solution  with  the  input  grid  generated  by  TCGRID  is 
very  stable. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A .   EXPERIMENT 

A  working  test  model  of  the  first  stage  of  the  SSME 
HPFTP  ATD  turbine  has  been  put  into  operation  in  the  Naval 
Postgraduate  School's  Turbine  Test  Rig  (TTR) .   The  TTR 
subsystems  have  been  successfully  exercised;  however,  all 
lines  and  connections  should  be  inspected,  cleaned  and 
repaired  as  necessary  before  any  measurements  are  taken  for 
performance  evaluation.   Connections  should  be  completed  to 
the  water  temperature  probes  and  torque  meter  on  the 
dynamometer  and  data  acquisition  should  be  automated.   Three 
methods  of  calculating  the  power  out  of  the  turbine  should 
be  automated  and  output  on-line:  namely,  using  the  flow  rate 
and  stagnation  temperature  drop  of  the  air  through  the 
turbine,  using  flow  rate  and  temperature  rise  of  the  water 
in  the  dynamometer,  and  using  the  rotational  speed  and 
torque  generated  in  the  dynamometer .   Once  the  system 
alignment  is  verified,  the  shroud  inserts  should  be 
installed  and  the  performance  measured. 

The  inlej  flow  circumferential  uniformity  should  be 
checked.   If  the  flow  is  not  uniform,  a  flow  straightener 
should  be  designed  and  inserted  into  the  inlet  duct.   It  is 
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also  recommended  that  the  inlet  and  exit  velocity  profiles 
be  determined  experimentally  . 

Ultimately,  the  turbine  outer  casing  will  need  to  be 
modified  to  allow  LDV  measurements  of  the  flow,  particularly 
in  the  tip  clearance  region. 

B.   NUMERICAL  SIMULATION 

A  working  CFD  model  of  the  SSME  HPFTP  ATD  turbine  first 
stage  stator  has  been  generated.   Four  solutions  were 
generated  on  two  different  computers.   It  is  recommended 
that  more  work  be  done  to  investigate  the  inlet -velocity 
profile  dependence  on  the  computer  system,  as  this  is  not 
acceptable.   An  investigation  is  needed  to  determine  whether 
the  anomaly  at  the  grid  cut  on  the  hub  is  grid  dependent,  or 
is  a  real,  physical  phenomenon. 

In  order  to  yield  more  realistic  flow  solutions,  the 
inlet  struts  that  hold  the  bullet  nose  upstream  of  the 
first-stage  stator  should  be  modeled  so  that  their  effect 
can  be  included  in  the  flow  solution.   Computations  can  be 
made  with  input  of  the  actual  inlet  velocity  profile,  once 
this  has  been  determined  experimentally.   For  an  analysis  of 
the  complete  stage,  the  rotor  geometry  must  yet  be  modeled. 
Such  a  modeling  of  the  downstream  portions  of  the  turbine 
will  allow  comparison  with  experimental  results  when  the 
turbine  is  tested  with  the  rotor  clearance  at  the  design 
value . 
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APPENDIX  A 


ENGINEER  T  NG  DRAW  I NGS 


Figure  Al .  Bearing  Outer  Spacer 
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Figure  A2 .  Probe  Mount 
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APPENDIX  B.   CALCULATIONS 

A.   MASS  FLOW  RATE 

The  following  assumes  the  9.0  in.  orifice  is  installed 
and  is  derived  from  Vavra  [Ref .  16] . 
First  solve: 


Now  find; 


where, 


m'^caY^, 


'tt 


\      Ti 


with  {=1; 


/f?/=cair1, 


~pJk 


N    ri 


*=0.2525  — 
z 


Z=1.9+0.24 ( 


ta 


100 


1) 


C=f(x) 


m=^m' 


in  lbm/s 


t1=  temperature  ahead  of  orifice  (°F) 

1^=  pressure  drop  across  orifice  (in.  H20) 

P1=    absolute  pressure  ahead  of  orifice  (in.  Hg) 
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1 .   Flange  Taps 


m^cccY*. 
\ 

PiK 

*i 

where  T1  =  t1  +  460  . 

Therefore, 

./        r          (t,-60)  ,  ,        hw, 
Jhf=6. 6445  [1+0.  0015 jL— — ]  (1-0.027-^; 


100 


~pJk 


Px    \     ^  +  460 


z=1.9+0.24  — -1) 

100 


x=0.2525  — 
z 


C=i+ 


0.0047 


77?=C-fn/. 


2 .   Vena  Contracta  Taps 


m 


^„ 


'=6  .7720  [1+0.0015 - ]  (1-0.  027  —  h 

100  Px  '\1    ti  +  460 


z=1.9+0.24  ( — —-1) 
100 


m' 

x=0.2525  — 

z 


C=i+ 


0.0029 


ih=^ih/ . 
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B.   TURBINE  POWER 


Power=mC  AT 
p 


Cp    =    0.24    Btu/lbm.'R    for    air 


AT=TT3-TT4 

TT4  =  Stagnation  Temperature  of  Air  at  Turbine  Exit  (*F) 
TT3  =  Stagnation  Temperature  of  Air  at  Turbine  Entrance  (*F) 
Since   1  Hp°s  =  550  ft.lbf   and   778  ft.lbf  =  1  BTU, 

Power=rh-  0  .  24  (  TT3-TT4)  221  . 
v         '  550 

C.   TURBINE  EFFICIENCY 

TT-  -TT 
Efficiency^ 


TT,-TTAl 


where  TT41  =  Isentropic  stagnation  temperature  is  given  by 

PT     Hi 
TT.,=TT,( 1)  y   . 
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APPENDIX  C.   COMPLETE  TABLE  OF  MEASUREMENTS 
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APPENDIX  D.   FORTRAN  PROGRAM  AND  INPUT  FILES 
RVC3D  Input  file 


'SSME  HPFTP  --  Indigo  --  Real  Slug  Flow' 

&nll  im=151  jm=31  km=45  itl=20  iil=70  &end 

&nl2  cfl=4.0  aviscl=0.0  avisc2=0.5  avisc4=1.5  ivdt=l  nstg=4 

itmax=985 

irs=l  epi=0.50  epj=0.60  epk=0.60  &end 

&nl3  ibcin=3  ibcex=3  isymt=0  ires=l  icrnt=50 

iresti=l  iresto=l  ibcpw=0  iqin=0  &end 

&nl4  emxx=0.212  emty=0.0  emrz=0.0  expt=0.0  prat=0.6500 
ga=l  .  4 

om=0. 000000  igeom=l  alex=-60.0  &end 

&nl5  ilt=2  tw=1.00  renr=6.651e6  prnr=.7  prtr=.9 
vispwr = .666  666 

srtip=0.0  cmutm=14 .  jedge=15  kedge=ll  iltin=0  dblh=0.00 
dblt=0.00  &end 

&nl6  iol=l  io2=151  oar=0.  ixjb=0  njo=l  nko=3 

jo=l  ko=13  17  21  &end 
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Inlet. f  -   Reads  RVC3D  output  and  builds  files  for  plotting 
Inlet  Velocity  Profiles. 

q*  *********************************************************  * 

c  inlet. f  reads  rvc3d  files  &  writes  ascii  files  for 

plotxy 

c     unit  1  =  input  xyz  file  (fort.l) 

c      unit  3  =  input  q  file    (fort. 3) 

c     unit  17  =  output  inlet  velocity  profiles 

c      unit  4  =  output  residual  history 

Q*    *********************************************************    * 

parameter (ni=152 ,nj=32, nk=46 ) 
real  x (ni , nj , nk) ,y (ni , nj , nk) ,  z (ni , nj , nk) 
real  qq (5 , ni , nj , nk) , resd (5000 , 5) 
real  vel (5)  , zz  (5) 
dimension  ii (5) 
c      i- values  are  hard-wired  below 
data  ii 76  0,70,75,80,90/ 

q*  **********   k  **********************************************  * 

c      read  q\ id  coordinates 

Q***********T.    *********************************************** 

read (1 , * ) im, jm, km 

read(l , *) ( ( (x(i,j,k) ,i=l,im) ,j=l,jm) ,k=l,km) , 

1  (((y(i,j,k),i=l,im),j=l,jm),k=l,km)/ 

2  (((z(i,j,k),i=l,im),j=l,jm),k=l,km) 

Q*   *********************************************************   * 

c      read  restart  file 

************************************************************ 

r  ead ( 3 , * ) imax , j  max , kmax 

read (3 , * ) f smach, alpha, re, time 
c 

icheck=iabs (im- imax) +iabs ( jm- jmax) +iabs (km -kmax) 

if ( icheck . ne . 0 ) then 

write(6,610) im,jm,km, imax , j  max , kmax 

stop 

endif 
c 

read (3 ,  *)   ((((qq(l,i,j,k)  ,i  =  l,im)  ,j=l,jm)  ,k=l,km)  ,1  =  1,5) 
c 
c      additional  residual  data 

r  ead (3,*)itl,iil, phdeg , ga , om , nr  e  s , dum , dum , dum , dum 

read (3,*) ( (resd (nr , 1) , nr=l , nres) ,1=1,5) 

Q*    *********************************************************    * 

c     vel  to  unit  17 

Q*    *********************************************************   * 

itr=im+l-itl 
ggm=ga* (ga- 1 . ) 
j=  jm 
c      normalize  x  by  chord 
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c      do  7  1-1, 5 

c     k=kk(l) 

c     xmin=x (im/2 , j , k) 

c      xmax=x(itl,  j,k) 

c     do  5  i  =  itl,  iti 

c     xmin=min (xmin, x (i , j , k) ) 

c    5  xmax=max (xmax, x (i , j , k) ) 

c     chord=xmax-xmin 

c     do  7  i=itl, itr 

c    7  x (i , j , k) = (x (i , j , k) -xmin) /chord 

c 

1  =  75 

j=jm 

k=km/2 

velinf =sqrt ( (qq (2 , i , j , k) **2+qq (3,i,j,k)**2 
l+qq(4,i, j,k)**2)/qq(l,i, j ,k)  ) 
c 
c     write (7 ,  305) (kk(l) , 1=1, 5) 

do  2  0  k=l,km 

do  10  1=1,5 

i=ii(l) 

veld)  =sqrt  (  (qq(2,i,  j  ,k)  **2+qq(3,  i,  j,k)**2 
l+qq(4,i, j ,k)**2)/qq(l,i,j,k) ) 
10  zz(l) =  z(i, j ,k) 
20  write ( 17 , 3  00) i,  (zz  (1)  ,vel (1) /velinf, 1  =  1, 5) 

300  format (i5, lOf 10 . 5) 

305  formate    k=  '  ,  5  ( 17x  ,  i3  )  ) 

610  f ormat ( '  *****    warning  *****•,/, 

1  '  im,  jm,  km,  read  from  input', 3i5,'  do  not 
match' , / , 

2  '  im,  jm,  km,  read  from  restart  file',3i5) 
stop 

end 
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program  shrink99 

This  inputs  data  and  multiplies  it  by 

:  .  1 1  I  ink   I  •  i  <  '  I  <  >  I  . 


0  .99  for  a  99% 


real  x (988)  ,p(988) 

open  (unit=2,file='hotlv.dat' ,status='old' ) 

open  (unit=15/file='lv.tcgin' , status = ' new1 ) 

read(2, *  )  (x  (i)  , i=l,988) 
do  20  i=l,988 

p  (  i  )  =  x  (  i  )  *  0  .  9  9 
20  continue 


the  data  out  in  blocks 
i-1, 88) 
i  =  89, 176) 
i=177 , 264 
i=265, 352 
i=353,405 
i  =  406 , 458 
i=459, 511 
i=512, 564 
i=565,617 
i=618,670 
1=671,723 
i=724  ,776 
i=777 , 829 
i=830, 882 
i=883,935 
i=936,988 
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TCGRID    Input    File 


>00(*OOO0O0C 


)  O  O  C?  O  O  r 


>  O  C>  O  O  ( 


???? 


"-????- 


I      l      I      I      i      I 


>  o  o  ( 1  r >  o  O  ' 


'???- -"  "• 


>  O  O  O  O  O  O  o  «  o  o  < 


»  O  o  o  r»  o  « 


>  O  O  o  © 


???"■"" 


■  o  o  — •  (n  a>  t 


o  o  ^  <"><"•  ' 


>  o  c  o  < 
jo  f-  «rt  r 
)  co        o  r-  v 


o  o  o  o  r 


& 


TV 


nor-  o  O  'J'  £•. 

O  "l  Ifl                —  rj  <t-  1-1  * 

0'><>  «  «>  »->  o  ' 

~-       .       .  v    hj  O  '-  •*  « 

«   0  9  ■"*        -  •  T.  o  ^  *•  « 
.   u   «,                                        ... 

«  j  o  'i  w  77*"" 


-ir-«BOt/>r-*^< 


>  «-*  o  O  f*  <*■  i 


P    -•    V     C    O  ct>     ■  1  f    T   »   r<   O  O  ^    "  "  O 

--■        «  0  /}  •  «j  **■  ~»  -»■  r-  n»  «•»     .     »    . 

— t-»         O  <T    ■  •  '*  <"»  O  O 


■  •»  o  c  t          •  >-  to 

Ti  ■-•  •    ■  fi  o  ▼    u  vfir-^9*c. 

.  <    I     tl  t*  L    ('  r     «-  .-"(rotr-f 

>.    1       :    t  »    ■  '  Hi-fronrrirOHNr-' 

e7!  fjE  .,-li-H-n 

•3    4  4    1          >  O  OT*>«^  On 
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EL       " 
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i  (O  O  (1  r<  '«  *r 


)  '3  O  O  C3  (  ID  C)  (J  O  O  '1  '  »  O  O  u 
>  ,oc*   «e  \r  .  (  n  o  m  "'  v  io  — *  '>«  i 

i  Of     <o  ■v  T'^ao'D  <*  hint  u'  -'oil1 


fl/llfl'fnOr*lft»H.I.1HHHH"4rtHi'IHHltHr'.4H.lf4rtN(«rl^     Mr-     .■■    ^    ■-*    «     ^     *f7     T 

nnon„01f^....,.,...—,.».«.," „„.... 
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•  •-•  O  O  O  ( 


)s»  O  ff)   "•  ■»  N  ' 
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